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In this work the shear viscosity of strongly interacting matter is calculated within a two-flavor 
Nambu-Jona-Lasinio model as a function of temperature and chemical potential. The general Kubo 
formula is applied, incorporating the full Dirac structure of the thermal quark spectral function 
and avoiding commonly used on-shell approximations. Mesonic fluctuations contributing via Fock 
diagrams provide the dominant dissipative processes. The resulting ratio rj/s (shear viscosity over 
entropy density) decreases with temperature and chemical potential. Interpolating between our NJL 
results at low temperatures and hard-thermal-loop results at high temperatures a minimum slightly 
above the AdS/CFT benchmark r]/s = 1/4-7T is obtained. 
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I. INTRODUCTION 


The quark-gluon plasma produced in heavy-ion col¬ 
lisions at RHIC and LHC is a hot and dense state of 
strongly correlated matter. It behaves like an almost- 
perfect fluid featuring a small ratio r\/s of shear viscos¬ 
ity to entropy density m- In this work we calculate 
viscous effects of interacting quarks within a two-flavor 
Nambu-Jona-Lasinio (NJL) model fJT-ITT]. A large-7V c 
scaling of the four-fermion vertex as inferred from QCD 
introduces a bookkeeping in which mesonic fluctuations 
(meson clouds around quarks) provide the dominant dis¬ 
sipative processes. The shear viscosity 77 is calculated 
using the Kubo formalism [T2] similar as in Refs. [13f 
[20] , The new element of the present work is that the 
full Dirac structure of the thermal quark self-energy is 
included when evaluating the Kubo formula, thus avoid¬ 
ing commonly used approximations p~7) I2TH23], At the 
same time we extend our previous studies in Ref. [2~4] . 

The present paper is organized as follows: in Section |TT| 
we discuss the NJL model from the perspective of large- 
N c scaling together with the gap equation and the Bethe- 
Salpeter equation. In addition we introduce the approxi¬ 
mation scheme for the quark-meson coupling used in this 


work. In Section III we develop the Kubo formula in¬ 
corporating the full Dirac structure of the thermal quark 
self-energy. In Section [TV| results for the quark self-energy 
Tj /3 generated by mesonic fluctuations are presented, as¬ 
suming first an on-shell approximation. In this case the 
coupling between quarks and mesons is dissipative only 
at sufficiently high temperatures where pions can decay 
on-shell into quark-antiquark pairs and thus the quark 
self-energy receives an imaginary part. In the next step 
constituent quarks off their mass shell are considered al¬ 
lowing for additional kinematic possibilities. Details of 
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the calculation are displayed in the Appendices [B] and [C| 
The results for the shear viscosity 77 and the ratio 77 /s 
are presented and discussed in Section |V} In the high- 
temperature region gluonic degrees of freedom become 
dominant and results from hard-thermal loop calcula¬ 
tions extend our NJL-model results [20]. Finally, Sec¬ 
tion [Vi] gives a summary of our most important findings. 


II. NJL MODEL AND MESONIC CLOUDS 

In this work the simplest two-flavor NJL model is used, 
including scalar and pseudoscalar interactions only: 

C = $ ( 1 $ - mo) ip + j [(t/u/O 2 + O/dTsT'f/’) 2 ] , (1) 

where = ( u,d) T is the isospin-doublet quark field, 
m 0 = diag(m u , nid) is the current-quark mass matrix 
(we work in the isospin limit, m u = = m 0 ), and t 

denotes the vector of three isospin Pauli matrices. The 
effective four-fermion coupling G is supposed to include 
non-perturbative (gluonic) dynamics. The large- N c scal¬ 
ing of QCD implies G ~ 1/N C in the NJL model. One 
should note that in QCD the color gauge symmetry is 
local whereas in the NJL model it is reduced to a global 
symmetry. 

In this large-lV c counting a hierarchy of Dyson- 
Schwinger equations can be introduced |251 125] , where 
the leading order 0(N®) is just the NJL gap equation in 
Hartree approximation: 

. _+ 

( 2 ) 

Lines with a black square denote full quark propagators 
whereas lines without denote bare quark propagators in¬ 
cluding the current-quark mass. In order to distinguish 
Hartree and Fock contributions generated by the four- 
quark vertex proportional to G we have introduced the 
wavy line in the last diagram. It indicates that a color 
trace is involved in the quark loop. Hence, the diagram- 
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Figure 1. Thermal masses of quarks (dashed curve) and 
mesons (solid curves) in the two-flavor NJL model with the 
parameter set given in Table [I] 


matic gap equation |2| reads 

m = mo — G^ip) . ( 3 ) 

It includes contributions from explicit chiral symmetry 
breaking, mo, and from the quark condensate: 

W> =- Y / d P t 1 “ U F ( E p) ~ U F ( E v)\ > 

7T Jo H/p 

(4) 

where n,p(E p ) = [1 + exp ((3(E p =F /r))] denotes the 
Fermi distribution functions with E p = \/p 2 + m 2 , and 
fi = 1/T is the inverse temperature. In the com¬ 
monly used mean-field approximation, the (thermal) 
constituent-quark mass is determined by solving just this 
Hartree part of the gap equation. The resulting quark 
mass within this approximation is shown in Fig. [T] 

At next-to-leading order the mesonic modes are ob¬ 
tained from the well-known Bethe-Salpeter equation 
(BSE) in random-phase approximation: 

=x + 

(5) 

When going beyond the mean-field approximation, the 
NJL gap equation includes Fock terms that are sup¬ 
pressed by 1 /N c : 


Input 

Output [MeV] 

mo 

G 

A 

m 

m-K 

U 

{ipi>) 1/3 

5.50 MeV 

10.1 GeV" 2 

651 MeV 

325 

140 

92.4 

-316 


Table I. NJL parameter set and resulting physical quantities 

The last (Fock) diagram in the extended gap equation 
© represents the mesonic clouds which couple mesonic 
fluctuations to constituent quarks. This coupling is 
described by Yukawa interactions with a single quark- 
meson coupling, g-Trqq, implied by chiral symmetry: 

AC = -g„ qq (jji'y 5 Tip • 7r + ipaip) . (7) 

Solving the Bethe-Salpeter equation ([5]) gives (renormal¬ 
ized) meson propagators from which meson masses can 
be extracted. The resummation of quark-antiquark scat¬ 
tering modes leads to meson propagators for any of the 
pions (P) or the sigma boson (S): 

D m = G+ GIi s / p D M = i _ ^ g/pG . (8) 

Here we have introduced the polarization tensors n s / p 
corresponding to the scalar or pseudoscalar quark- 
antiquark loops. They have the form 

H S/,p (p, u> n ) = 8 N c h + 4 N c N s / p I 2 (p,co n ) , (9) 

where N p = — (co 2 + p 2 ) and N s = N p — 4m 2 refer to 
the pion and sigma modes, respectively. Explicit expres¬ 
sions for I\ and l 2 (p,uj n ) can be found in Appendix [A] 
Note that the gap equation in Hartree approximation pi 
involves also the loop integral Ip. 

m = mo + 8GN c mI± . (10) 

Poles of the meson propagators Pm can appear only in 
Minkowski space, therefore one performs the analytical 
continuation iw„ i —> oj + ie. In general the polarization 
tensor is complex which is evident at high temperatures 
where a mesonic resonance instead of a bound state is 
realized. Therefore, we define the meson mass tom at 
p = 0 as a solution oQ 

Re (Oj — mtim) = 0 • (11) 





( 6 ) 


A self-consistent treatment of the gap equation including 
Fock terms with mesonic modes and the Bethe-Salpeter 
equation describing these mesonic modes is approximated 
by a common procedure mi nz]: the mesonic fluctuations 
are evaluated using the Hartree solution of the gap equa¬ 
tion only which ensures a consistent large--/V c counting of 
the NJL model up to next-to-leading order 1/N C . 


In the pseudoscalar channel the (renormalized) pion 
propagator reads 


E TZ ( p , CJ n ) 


G 

^+4GN c (co 2 n +p 2 )I 2 (p,u Jn ) ’ 


( 12 ) 


1 In principle, there is a momentum dependence of the pion mass. 

We have chosen to define masses always in the reference frame 
of the heat bath, so p = 0 and 
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Figure 2. Momentum dependence of the pion mass at two 
different temperatures 



from which the pion mass can be calculated by solving 
mo 1 


ml = 


m 4 GN C Re J2(0, —im n ) 


(13) 


If one considers the scalar channel instead, the mass of 
the sigma boson can be extractec0 

m 2 = m 2 + Am 2 . (14) 


At high temperatures T, far above the chiral transition 
temperature, the scalar and pseudoscalar modes tend to 
degenerate: m 2 m 2 . This goes along with the restora¬ 
tion of chiral symmetry. The results for the thermal 
meson masses are also shown in Fig. |T] The parame¬ 
ter set which has been used to calculate m(T), ?n w (T) 
and m a (T ) is given in Table [I] It reproduces physical 
values for m ^ = 140 MeV, for the pion decay constant 
= 92.4 MeV, and realistic values for the constituent- 
quark mass m = 325 MeV and the chiral condensate 
('tpip) = —(316 MeV) 3 . Due to the absence of confine¬ 
ment in the NJL model, the thermal pion can decay on- 
shell into two thermal constituent quarks. The critical 
temperature for this is called Mott temperature 7m and 
determined by ?tv(7m) = 2to(7m). At vanishing chemi¬ 
cal potential it has the value 7 m ~ 212 MeV. Therefore, 
for the shear viscosity arising from mesonic fluctuations, 
on-shell dissipative effects are possible only for T > 7m. 

We continue with examining the validity of the stan¬ 
dard pole approximation for determining meson masses. 
In Minkowski space, the solution of the BSE (at finite T 
and p) reads for the pionic mode: 


Dtt(p , -iw) 


G 

— 4GIV C (w 2 — p 2 ) I-2( P , —iw) • 


(15) 


2 Note that this relation between and m n is valid only if the 

energy dependence of I 2 is negligible. For high T pions and the 
sigma boson are degenerate and one has 72 (—i?7v) ~ / 2 ( — ima). 


The fact that hip, — icd) is energy and momentum depen¬ 
dent implies that the standard pole-mass approximation, 


^Vqq, static 

w 2 — p 2 — m 2 + is ’ 


(16) 


does not reproduce the full (p, oj) dependence of the prop¬ 
agator. A general quark-pion coupling can be introduced 

by 


5 ? qq (w,p) = -(w 2 -p 2 -m 3 (p))T7 7r (p,-iu;) , (17) 

with the momentum dependent pion mass (p) defined 
as a solution of 


ReD^ x (p, —im w (p)) = 0 . (18) 


The results for the momentum-dependent pion mass 
m^ip) are shown in Fig. [2] The pion becomes more mas¬ 
sive when it carries additional momentum. This quali¬ 
tative behavior of m^^p) is consistent with the fact that 
the constituent-quark mass rri(p) decreases as function 
of momentum. The associated tendency towards chiral 
symmetry restoration weakens the Goldstone boson char¬ 
acter of the pion. 

As usual we define the quark-meson coupling as the 
residue of the full meson propagator at vanishing mo¬ 
mentum |28| : 


^qq = - ^2 ^ (°, 


(19) 


From the pion propagator D in Eq. (151 we get imme¬ 
diately 


Su-qq =4Vc 7 2 (—iw) + U! 


duj 2 ) 


= 9 


-2 


1 + 


d In 


Trqq,static I ^ / 2 (_kj) dw 2 


( 20 ) 


We have identified the static quark-meson coupling, 
f? 7 rqq,static? where the energy dependence of J 2 (0,-iw) 
is neglected and find 3~ q 2 qjStatic = 4V c / 2 (0, -im n ), cf. 
Eq. (16). In Fig. [3] we compare the two approaches for 
calculating the quark-pion coupling. When staying in a 
50% interval around the pion pole, 0.5771^ < ui < 1.5 m n , 
we find indeed that the usual treatment, fixing g^qq at 
the pion pole, is a rather good approximation. It is in¬ 
teresting to note that for vanishing three-momentum at 
T = 150 MeV the approximated quark-pion coupling, 
<Aqq,static overestimates the actual coupling when evalu¬ 
ating at energies smaller than the pion mass. In contrast, 
for high momenta but the same temperature, the cou¬ 
pling is underestimated by up to 6% which is still a good 
approximation. Only at very high temperatures and far 
away from the actual pole mass sizable deviations oc¬ 
cur. However, apart from such extreme values, one can 
conclude that the approximated quark-pion coupling is 
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Figure 3. Corrections to the quark-pion coupling in a 50% 
range around the pole position, (a) Vanishing pion momen¬ 
tum; (b) pion momentum p = p = 400 MeV 


acceptable and corrections beyond Eq. (20 1 contribute 
at the order of only a few percent. Therefore, we do 
not take any momentum dependence into account apart 
from the derivative correction in Eq. (20 1 . This leads to 
a simplified treatment of the mesonic fluctuations since 
the quark-meson coupling is just a constant and therefore 
does not affect the momentum integration. 


III. KUBO FORMALISM AND QUARK 
SELF-ENERGY 

In the Kubo formalism the shear viscosity is related to 
a correlator of the energy-momentum tensor. Assuming 
an infinite homogeneous medium close to thermal equilib¬ 
rium, the (frequency dependent) shear viscosity is given 

by M- 


v{u) = j d te' ut J d 3 x (T /JjV (t,x),T fJ,J ( 0,0)), 

( 21 ) 

with the energy-momentum tensor of quarks: 

d £ 

T ix V = ^ ~ = tylp&A ~ ■ ( 22 ) 


The correlator in the integrand of the Kubo formula (21 1 
is defined through the thermal expectation values with 

(A, Y) = (Ae-« ff Fe«% , (23) 

with H the Hamiltonian. As discussed in our previous 
work |24| . the shear viscosity can be written in terms of 
the quark spectral function p = — ^ Im Gr (with Gr the 
retarded quark propagator, see also Refs. |211l^5] ): 


vL=o = 




T 


x Tr [72 p{s,p) 72 p{e,p)\ ■ 


(24) 


In Ref. |2U the shear viscosity i?[r(p)] has been ex¬ 
plored assuming a simple parameterization of the thermal 
quark propagator with a schematic (momentum depen¬ 
dent) spectral width, T(p;T, p): 


Gr(po,p) = 


1 


ft — m + i sgn(p 0 )T(p) 


(25) 


from which the quark spectral function can be derived. 
The general Dirac structure is, however, richer than the 
parameterization in Eq. (25 ). Due to the breaking of 
Lorentz invariance in the thermal medium, three func¬ 
tions A 1 B , G are necessary to specify the quark spectral 
function: 

p(po,p) = [mA + po^oB -p- 7 C], (26) 

with a denominator D. These four functions depend on 
the (off-shell) energy po, the three-momentum p , and the 
thermal parameters T and p. They can be determined 
from the thermal quark propagator in Minkowski space 


Gr(po,p) = 


1 


i — TO — £ 


m(l + So) + Po7o(1 + £ 4 ) ~ P ■ 7(1 + £ 3 ) 
P§( 1 + S 4 ) 2 - P 2 (l + S 3) 2 - m 2 ( 1 + £ 0 ) 2 


(27) 


The general Dirac structure of the thermal quark self¬ 
energy is 


£ = wE 0 + p ■ 7 E 3 - P 070 S 4 , 


(28) 


with three dimensionless functions Y,j ( po , p ). Incorpo¬ 
rating mesonic fluctuations within the NJL model, the 
self-energy £ receives contributions from three pions and 
one sigma boson: 


£, = 3EJ 


S-, for j = 0,3,4 , 


(29) 


with E)’ //I> further specified after Eq. (38). 

For the present calculation we take into account only 
the relevant imaginary parts, 


Im E j = Pj . 


(30) 






















5 


In doing so, we ignore the momentum dependence of the 
constituent-quark mass as it arises from mesonic Fock 
contributions. Formally, this approximation is equiva¬ 
lent to readjusting the NJL parameters and introducing 
a new set (too, G, A) new that will depend on the thermal 
variables T and p, and on energy and momentum. The 
resulting thermal quark propagator is: 


Gr(po,p) 


TO (1 + ipo) + Po7o(l + 1 Pa) 
N\ + 2i N 2 


P ' 7(1 + 1 P 3 ) 
(31) 


with the two auxiliary functions Ni^{po,p) : 

Ni = pl{ 1 - pi) -P 2 { 1 - pi) - rn 2 ( 1 - pi) , 

N 2 = plp4 - P 2 P3 - rn 2 p 0 . 


(32) 


The four energy and momentum-dependent functions 
A,B,C,D parameterizing the quark spectral function p 
in Eq. (261 are thus identified as: 


A — poNi — 2 N 2 , B — P4N1 — 2N2 , 
C = P3 N x - 2 N 2 , D = N 2 + AN 2 . 


(33) 


The evaluation of the shear viscosity is now reduced to 
carrying out all the traces for the integrand of Eq. (24): 

AN c N f 


Tr [72/572P] = 


2 D 2 


[-m 2 A 2 +p 2 0 B 2 -p 2 C 2 + 2 p 2 y C 2 ] . 

(34) 


After angular integration the shear viscosity reads: 


V = 


2 N C N { 

3tt 3 T 


/ oo nA 

de dp rip (e)(l — rip (e)) 

-00 J 0 


D 2 {e,p) 


-m 2 A 2 (e,p)-yC 2 (e,p) + e 2 B 2 (e,p) 
5 


(35) 

with A = 651 MeV the NJL cutoff and N C N{ = 6 . 
Quite remarkably, negative and positive contributions 
balance to an overall positive shear viscosity 77 > 0. It 
is important to note that according to the representa¬ 
tion in Eq. (35), 77 is an even function of the chemi¬ 
cal potential p. This is ensured by a separate integra¬ 
tion over positive and negative energies and the property 
I(—e,p) = /(e, —p) of the entire integrand. 


IV. QUARK SELF-ENERGY FROM MESONIC 
FLUCTUATIONS 

In this section we evaluate the quark self-energy aris¬ 
ing from the Fock diagram with mesonic fluctuations 
([6]). They introduce non-vanishing imaginary parts pj 
at next-to-leading order in l/N c . 


A. On-shell quarks 

Let us first consider quarks with on-shell kinematics 
Po = m 2 +p 2 ■ The Matsubara frequencies for in-medium 


1.0 

0.5 

0.0 
w" 

6 -°.5 

7 

- 1.0 
-1.5 
- 2.0 

0 50 100 150 200 250 300 

p [MeV] 

3 
2 
1 

W 0 


-2 
-3 
-4 

0 50 100 150 200 250 300 

p [MeV] 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

T = 260 MeV, p = 0 ( 

m(T, p) = 35 MeV 

- 

- -Im I 0 


- -Im I, 

/ 

. -Im I 4 

/ 

;r - case : 


: 

T = 220 MeV, p = 0 
m(T , p) = 72 MeV : 



__ „ r " Tr ‘ ” 


-Im E 0 
-Im Z 3 
-Im E 4 


n — case 


Figure 4. The (negative) imaginary parts of the quark self¬ 
energy contributions E j, j = 0,3,4, from mesonic fluctua¬ 
tions. They have been defined in Eq. (381. 


quarks are v n = (2 n + 1 ) 7 rT — ip. Note that the frequen¬ 
cies for antiquarks are u* = ( 2n + l)nT 4- ip. There are 
N 2 — 1 = 3 equal contributions from the pseudoscalar 
channel (pions, F p = iys) and one contribution from the 
scalar channel (sigma boson, T s = 1). The correspond¬ 
ing self-energies are calculated as: 

yf P (p, i'n) = * » 

= (36) 

x r s/p G p (p - q,v n - v m ) , 
with the thermal quark and meson propagators, 


G}(p,v n ) 


Gp{p,u} n ) 


Vnli -p-~Y + m 
l' 2 +p 2 +m 2 
1 

U} 2 + p 2 + TO . 2 


(37) 


where tu n = 2 mrT are bosonic Matsubara frequencies. 
Its Dirac structure has the following form: 

S |g /P (p, Vn) = E 0 - p ■ 7 s 3 + t /„74 E 4 , (38) 
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with three dimensionless functions Yjj(p,v n ). The plus 
and minus sign in front of Eo refers to the (scalar) sigma 

boson and the (pseudoscalar) pion, respectively: E 3 ' 4 = 
S/P 

E 3 4 but E 0 = ±Eq. We note that in the single-width 


approximation in Eq. (25) would give ImEo = — ^T(p) 
and E 3 4 = 0 . 

We now analytically continue the quark self-energy to 
Minkowski space, v n > —ipo, and extract the imaginary 
parts of Ej relevant for calculating the shear viscosity. 
The detailed derivation can be found in Appendix [B] 
Here we only state the results (j = 0,3,4): 


Im Ej(p,-ip 0 ) 


?7rqq 


167Tp J E 

1 " J-'Tt 


d EfJFj [n B (E b )+n F (E f )] , 


(39) 


These results for Im Ej(p, — ipo), j = 0,3,4, will be 
used for the evaluation of the shear viscosity (35), where 
Po = \Jm 2 + p 2 is a function of the quark momentum 
due to the on-shell treatment. Their dependence on the 
momentum p is shown in Fig. [4] for the pion case at two 
different temperatures, T = 220, 260 MeV, and vanishing 
chemical potential. Due to the explicit analytical form 
of the self-energy contributions, they can be easily imple¬ 
mented and numerical issues arise only from handling the 
peaks in the integrand of the Kubo formula (351, cf. also 


Ref. |2'4] . We emphasize that its energy and momentum 
integrals are carried out independently and the functions 
A, B , C, D in its integrand remain off-shell. 


B. Off-shell quarks 


with E b = Ef + p 0 , po = \Jm 2 + p 2 , the weight factors 
E 0 = l, 

Tom — 2 m 2 — 2Efpo 


E 3 = 


Ea = — 


2 p 2 


(40) 


El 

Po 


and the integration boundaries 
1 


E, 


max, min 


2 m 2 - 


(tom - 2to 2 ) yjm 2 +p 2 
±p m M \jmu - 4?n 2 


(41) 


The index M is either S (sigma boson) or P (pion). The 
remaining integral over Ef in Eq. ( |39| ) can be performed 
and one finds the following analytical expressions: 

Im EoCp, -ipo) = ^51 Tin (jBmin) " b(£U “ +Po) 


167 rp n F (£ max ) n B {E min + p 0 ) ’ 

(42) 

with n F denoting the antiquark distribution function de¬ 
fined after Eq. @ and n B {E) = [exp((3E) — 1] 1 the 
Bose distribution function. 

The representations of E 3 and E 4 contain energy- 
dependent prefactors, J -3 and J- 4 , respectively, which lead 
to more complex results. After introducing the auxiliary 
function 


n(E) = (E + p 0 ) In n~(E) 

1 


-TLi 2 - 


n B (E + p 0 ) 


-T Li 2 1- 


n F (E) 


(43) 


one obtains: 


ImE 3 (p, -ip 0 ) = 


2 p 2 +m 


M 


2 p 2 


Im En 


9lmPv T n ^ E )\f 


q 2 T 

Im E 4 (p, -ip 0 ) = Im E 0 + U{E)\^ 


I 67 rp 3 


167TPPQ 


\E„ 


( 44 ) 


So far we have treated the external quark in the Fock 

S /P — 

self-energy Ef (p, — ipo) in Eq. (36 1 as an on-shell par¬ 


ticle with Po = m 2 + p 2 when determining the imaginary 
parts of Ej (p, —ip 0 ), j = 0, 3,4. According to the general 
Kubo formula for the shear viscosity derived in Eq. (241, 


the quark spectral function p(po,p) enters for off-shell 
kinematics. Using on-shell expressions for Im E j is a 
commonly used but unnecessary approximation. While 
in the on-shell approximation only one dissipative process 
(meson decay into a quark-antiquark pair) gives rise to 
an imaginary part, the off-shell situation features several 
dissipative processes. In Appendix 0 the detailed ana¬ 
lytical calculation of the off-shell imaginary parts of the 
quark self-energy from mesonic fluctuations is presented. 
Here we state only the results: 


Im Eo ff = 


Im E? ff = 


cf 7rqq 

167TP 

rqq 

167TPPQ 


(JI + JII + JIH) ; 

(. K 1 + K u + K 1U ) 


(45) 


and 


Im ES ff = 


m M +P 2 ~ Po ~ m 2 

V 


Im Eg ff + ^ Im E 4 

pz 


Po 


off 


(46) 

It is interesting to note that in contrast to the on-shell 
results in Eqs. (42) and ( [44] ), the off-shell imaginary parts 
feature a vacuum contribution included in J m and A ' 111 
which do not vanish in the limit T, p — > 0. 


V. RESULTS FOR THE SHEAR VISCOSITY 

We are now ready to present results for the shear vis¬ 
cosity 77 in the NJL model. First we use the on-shell 
expressions for the imaginary parts Im E j, j = 0,3,4, 
written in Eqs. (421 and (44) to evaluate the Kubo for¬ 


mula (35) numerically. The temperature dependence of 
the viscosity is shown in Fig. [5] for two values, /i = 






























7 



Figure 5. Temperature dependence of shear viscosity calcu¬ 
lated from the NJL model in its large- N c expansion for van¬ 
ishing quark chemical potential and p = 200 MeV. See the 
discussion in the text. 



Figure 6 . Temperature dependence of the ratio 77 /s for van¬ 
ishing quark chemical potential and p = 200 MeV. See the 
discussion in the text. 


0, 200 MeV, of the quark chemical potential. Due to the 
on-shell restrictions only the temperature range above 
the Mott temperature, T > Tm, is accessible. One has 
Tui/j, = 0) = 212 MeV and T u (p = 200 MeV) = 
171 MeV. We observe an overall decreasing function rj(T) 
and also decreasing values 77 ( 77 ) for increasing the chem¬ 
ical potential. A small shear viscosity reflects a highly 
correlated system: stronger interactions with the thermal 
medium lead to a lower value of 77 [23 . We conclude that 
the quark plasma described by the NJL model, where the 
shear viscosity is induced by mesonic fluctuations occur¬ 
ring at order l/N c , becomes more strongly correlated for 
both increasing temperature and chemical potential. 

Now we turn to the ratio 77 /s, shear viscosity to entropy 
density. Consistent with the l/7V c -approach, we use for s 
the entropy density of non-interacting constituent quarks 
with (T, /^-dependent masses: 

AT 7VT /»A,00 

s(T,fi)= c 2 f / d pp 2 [— In rip (E) — In rip (E) 

71 Jo 

+/3(E + p)n(E) + (3(E- p)n,p (E)] , 

(47) 

where the upper boundary of the momentum integral, 
/ A ’°°, encodes E = y/m 2 + p 2 for p < A but E = 
\Jm,Q + p 2 for p > A, with the current-quark mass mo 
instead of the constituent-quark mass m. This so-called 
soft-cutoff scheme ensures the correct Stefan-Boltzmann 
limit of s(T, p) at high temperatures. Inspection of Fig. 6] 
shows that the overall scale of the ratio 77 /s is comparable 


to 1 / 47 T [291130] . However, for large enough temperatures 
it undershoots the AdS/CFT benchmark as it can happen 
also in other quantum field theoretical models [314433] . 
In the NJL model this happens for vanishing chemical 
potential at Tw 275 MeV, and for a finite chemical po¬ 
tential, p = 200 MeV, at a somewhat lower temperature 
T « 260 MeV. 

Furthermore, we compare our results for 77 /s to those 
from lattice QCD, [M] ;35j, which are shown as squares 
with error bars in Fig. [ 6 ] They have been obtained 
within pure-gauge QCD and suggest a rising ratio 77 /s 
for T > 250 MeV, a behavior which is not found in the 
NJL model. This qualitative difference can be explained 
by considering results from hard thermal loop (HTL) cal¬ 
culation in QCD [SHECZl- At leading logarithmic order 
one finds the behavior [38] : 


C1T 3 

ocl In (C 2 /a s ) 


(48) 


with flavor-dependent coefficients C\ and C 2 . Conse¬ 
quently, the dimensionless ratio 77 /s scales as 77 /s ~ a ~ 2 
at leading order. For increasing temperature the QCD 
coupling becomes weak, a s —> 0 , and the ratio 77 /s rises 
with T according to the HTL results. Lattice QCD 
suggests that this trend sets in at rather low tempera¬ 
tures, where HTL calculations are not applicable since 
they are based on perturbative-QCD and resummation 
techniques. The main reason for the rising behavior 
of 77 /s in lattice QCD are the weaker correlations be¬ 
tween the gauge bosons towards asymptotic freedom. 





















Figure 7. Individual and summed ratios t]/s(T) from the NJL 
model (low-T region) and from HTL calculations (higli-T re¬ 
gion) at vanishing quark chemical potential. See the discus¬ 
sion in the text. 


In contrast to this, the NJL-model coupling G remains 
constant and the viscous effects from mesonic fluctua¬ 
tions are growing in the considered temperature range 
180 MeV < T < 300 MeV. As a consequence, the NJL 
model provides 77 and 77 /s decreasing with T and 7 1 . Note 
that for large T the thermal quark mass, m(T) ~ gT, 
dominates the constituent-quark mass derived within the 
NJL model. In our results the quark mass becomes small 
in this temperature region: m —> mo, cf. Fig. [T] 

The open circles in Fig.[ 6 ]are the results from Ref. |39j . 
where the shear viscosity has been evaluated from the 
basic Kubo formula (2l]) using cross sections utot from a 
parton cascade model with elastic two-body collisions for 
gluons only. Their results are described by 


77 _ 0.195 
s a totT 2 


(49) 


and utot = 9 mb = 0.9 frn 2 has been used to obtain the 
open circles in Fig. [ 6 ] In comparison to our NJL results 
one gets a decreasing but flatter ratio 77 /s in that ap¬ 
proach. The assumption of a temperature-independent 
total cross section a to t does not describe the high-T 
behavior of HTL calculations and suggested by lattice 
QCD. 

The rising behavior of 77 /s has been parameterized in 
Ref. [20] as 


HTL 


a. 


( 50 ) 


with a = 0.2 and 7 = 1.6 extracted from a combined fit 
to results from functional-renormalization-group meth¬ 
ods and HTL calculations. Note that 7 « 2 as expected 
from the pure HTL result for gauge theories, Eq. (48 1 . 
Ref. [20] has used the following form for the temperature 
dependence of the running QCD coupling j40fl42] : 


s(T) = 


47t z 2 — 1 
/3q z 2 lnz 2 ’ 


(51) 


with the beta-function coefficient 0q = (111V C — 27Vf)/3, 
and the reduced temperature z = 0.79 T/T c , where T c ss 
155 MeV. 

The shear viscosity 77 obtained from HTL calculations 
is induced by dissipative processes in the gauge sector, 
whereas 77 calculated from the NJL model arises from 
mesonic fluctuations in the quark sector. We interpolate 
between the low-T (NJL) and high-T (HTL) domains by 
taking the sum of the two corresponding ratios 77 /s, as it 
is suggested in Ref. | 2 U| . The resulting summed ratio is 
shown as the solid line in Fig. [7] It develops a minimum 
at T m ; n = 295 MeV with 77 /s(T m ; n ) = 0.29 > 3.6/47T due 
to the change between quarks and gluons as active de¬ 
grees of freedom. In comparison to the analogous results 
T^f = 200 MeV and 77 /s(T^°) = 0.17 from Ref. [2D], 
both the minimal value of 77 /s and its location are shifted 
to higher values in the present work. The main reason 
for this is the rather high chiral crossover temperature, 
T c = 190 MeV, in the two-flavor NJL model. Clearly, 
Figure [^should be taken just for qualitative orientation. 
The position and value of the minimum in 77 /s depend 
sensitively on the detailed interpolation between the low- 
T and high-T domains. Taking the sum of the corre¬ 
sponding ratios at low and high temperature is only one 
possible way. 

We have also numerically evaluated the Kubo formula 


for the shear viscosity (35) using results for the off- 


shell spectral functions pj = Im E j = 0,3,4 given 
in Eqs. ( |45| and ( [46] ). The corresponding results are 
shown in Fig. [ 8 ] and compared. The new feature of the 
proper off-shell treatment is the presence of viscous pro¬ 
cesses in the whole temperature region. Therefore, the 
constituent-quark mass does not provide any restriction 
on a finite shear viscosity. It is interesting to observe 
that at p = 0 the results at small T smoothly join those 
of the on-shell approximation where the Mott-condition 
Tom > 2 ?77 has to be fulfilled. The quantitative difference 
in this region is almost negligible. At p = 200 MeV the 
difference is an almost constant factor shifting the viscos¬ 
ity to higher values, but the overall behavior of rj(T, p) 
is not changed. 

One can explain this qualitative agreement by the 
peaking of the integrand in the Kubo formula: the main 
contribution to 77 in Eq. (35) is collected around the min¬ 
imum of the denominator D(e,p), cf. Fig. 2 in Ref. |24j . 
This essentially leads to the on-shell approximation as 
one can argue with Eqs. ([32]) and (33): 


D —»• 0 <t=> Ni , N 2 0 => Pq — p 2 — m 2 = 0 . 

(52) 
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T [MeV | 



Figure 8. Comparison between the on-shell and off-shell cal¬ 
culation of the shear viscosity at vanishing quark chemical 
potential (left panel) and p = 200 MeV (right panel) 


We conclude that the off-shell treatment provides only 
subleading corrections due to the peak structure of the 
integrand of the Kubo formula. One would expect that 
the shear viscosity becomes smaller in the off-shell treat¬ 
ment since more dissipative processes are at work. But on 
the contrary, the complicated arrangement of the imag¬ 
inary parts pj, j = 0,3,4, in the integrand of p leads 
eventually to an increasing shear viscosity compared to 
the on-shell approximation. 


VI. SUMMARY AND CONCLUSION 

In this work we have investigated the shear viscosity of 
hot and dense quark matter described by a large- N c NJL 
model for two flavors. We have used the Kubo formal¬ 
ism and calculated the shear viscosity p from a thermal 
quark spectral function with inclusion of its full Dirac 
structure. Instead of a single width there are now three 
(off-shell) imaginary parts which determine the positive- 
definite shear viscosity. 

In the large-_/V c counting, the dominant dissipative 
process arises from mesonic fluctuations. They are dy¬ 


namically generated by virtual quark-antiquark loops re¬ 
summed to all orders in the non-perturbative Bethe- 
Salpeter equation. The mesonic Fock contribution to the 
gap equation are of subleading order l/N c . We have cal¬ 
culated the three components £j of the quark self-energy 
provided by the mesonic Fock term, both for on-shell 
and off-shell kinematics. Evaluating the Kubo formula 
with the input pj = Im £j, we have found a decreas¬ 
ing shear viscosity as function of both temperature and 
quark chemical potential. At vanishing chemical poten¬ 
tial, the proper off-shell treatment extends the on-shell 
approximation into the low-temperature region where the 
on-shell viscous effects are kinematically excluded. Apart 
from this, off-shell corrections have no further quantita¬ 
tive or qualitative influence. However, at finite quark 
chemical potentials, off-shell effects shift the shear viscos¬ 
ity to higher values but its overall qualitative behavior is 
not changed. 

We have observed that the dimensionless ratio 77 /s un¬ 
dershoots the AdS /CFT benchmark 1/An at large enough 
temperatures in the NJL model. Combining our results 
for the shear viscosity with perturbative results from 
hard-thermal-loop calculations in the high-T region, we 
find that the ratio p/s develops a minimum well above 
the AdS/CFT benchmark. The interpolated results com¬ 
pare reasonably with those from lattice QCD regarding 
the overall behavior and scale of the ratio p/s. However, 
since the chiral crossover temperature in the two-flavor 
NJL model, T c ss 190 MeV, is larger than the lattice- 
QCD result, T c ss 155 MeV, the onset of the dominant 
viscous effects is shifted to higher temperatures by the 
Mott condition. 

In summary one can conclude that the correlated quark 
matter described by the NJL model features a small ratio 
shear viscosity over entropy density as it is characteristic 
for a perfect fluid. 
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Appendix A: Meson propagators from the BSE 

The meson masses are derived from the Bethe-Salpeter 
equation (BSE) as pole-masses of resummed quark- 
antiquark scattering modes: 

D M = G + GU s / p D M = i _ ^ s/pQ , (Al) 
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where II s / p denotes the polarization tensor: 

n s/p (p, Wn ) 

= 8 N C T Y 


d 3 q =F m 2 + v m (y m - u n ) + q{q - p) 


( 27r ) 3 Wm + E q][(^m - W„) 2 + E 2 a ] 


= 4 N C T Y 


d 3 q 


(2t r) 3 


1 


1 


+ 4 N C N S/P T Y 




"m + E q {Vm - U n ) 2 + E\ _ 

d 3 q 1 

(2tt)3 [u 2 m + E 2 }[{v m -u n ) 2 + El\ ’ 

(A2) 


with the energies E 2 = q 2 +m 2 and = ( q — p ) 2 + m 2 . 
It can be expressed as 


n s ^ P (p,oo n ) = 8 N c h + 2N c N s / p I 2 (p,u n ) , (A3) 


where we have denoted the momentum-independent part 
of the tensor by 1 i, whereas the momentum dependence 
is encoded in l 2 (p,co n ). In addition we have introduced 
N p = — (w 2 + p 2 ) and N s = N p — Am 2 describing the 
pion and sigma modes, respectively. We have defined: 


I 2 {p,oj n ) = TY 

d 3 q 


d 3 q 

( 27r ) 3 Win + E q][(^n - V m ) 2 + E 2 a \ 

2 E + {iui n - E + ) (n F (E q ) + n£(E A )) - (iw n + E + ) (n^(E q ) + rip (E A )) 


(27r) 3 4 E q E A 


“l + E \ 


“l + E l 


+ 


(iw n + E_) (ri£(E q ) - n£(E A )) - (iu n - £_) (n F (E q ) - n F (E A )) 


oo. 


2 + E 2 

rt 1 — 


(A4) 


having defined E± = E q ± E A . The momentum- Appendix B: Details of the on-shell calculation 

independent part of the polarization tensor reads 

The individual contributions are given by 


t = t s t [ d3g _ 1 

1 2 ^J ( 2 tt ) 3 vl + q 2 + to 2 

= 47T 2 J 0 d ^£T ^ ~~ n F^ E v) — n F ( E p)) ■ 


(AS) 


As mentioned already before, Ii refers also to the chiral 
condensate, which can be seen from inspecting Eqs. |3| 
and (flOj) : 


m — mo (V’V’) 

8 GN c m 8 N c m 


(A6) 


In the calculation of II s / p one has to take care of sev¬ 
eral minus signs: the fermion loop gives a global minus 
sign, in the pseudoscalar channel one has i 2 = — 1 and 
additionally { 75 , 7 ^} = 0 and { 7 *, 7 ^} = — 2 f )., 7 for the 
Euclidean gamma matrices. 


So — g^ qq T Yj 
S 3 = glqqT Y 
^4 = glqqT Yj 


d 3 q 


(27r) 3 + E 2 (v n - v m ) 2 + E 2 ’ 

d 3 q p q 1 1 

(27r) 3 p 2 + E 2 f (v n - u m ) 2 + E 2 

d 3 q v m 1 _1_ 

(27t) 3 v n + E 2 f (y n - v m ) 2 + El 

(Bl) 

with the energies E 2 = q 2 + m 2 and El = (p — q) 2 + 
As always, the Matsubara sums can be carried out 
leading to some finite result with a combination of Bose 
and Fermi distribution functions. We arrive at 




So,3 (p,Vn) = glqqj 

E~\~ Z\ E— Z 2 


d 3 q 

(2tt) 3 


Ti 


0,3 


1 


2 EbEf 

m, z 3 


E %+< E 2 -+<) ( E l+^)( E 2 -+<). 

(l’>2) 

with E± = Eb ± Ef. The quark-meson coupling, gMqq, 
can be pulled out of the integral since no momentum 
dependence is taken into account as it has been justified 
before. We have introduced Tq 3 as 


To = 1 


•F 3 = 


p q Tom - 2m 2 - 2Efp 0 ( B3 ) 

p 2 2 p 2 
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and have denoted the combinations of Bose and Fermi 
distributions as Zi(Eb, Ef): 


’(Ef)), 


z i = 1 +n B (E b ) - - (r4(Ef) - 
Z 2 = n B (E b ) + i (n^(E f ) + n^(E f )) > 0 , 

Z 3 = n F (Ef) — n F (Ef) > 0 . 

The Bose and Fermi distributions read 

n n( E ) = ef}E _ 1 > n-p(E) = n F (E^p) = e p( E ^») + 1 > 

(B5) 

where the ± signs denote quark and antiquark distri¬ 
bution functions, respectively. When carrying out the 
Matsubara sum also for the £4 part of the self-energy, 
we get: 


and converted the momentum integral to an energy inte¬ 
gral using qdq = Ef d Ef. The ill-conditioned 0 term can 
be removed by the following consideration: from Eq. (B91 
it is clear that |£| < 1 is fulfilled if and only if 


- 1 < 


rr 2 ^2 „2 

E b ~ m M ~P - 9 

2 pq 


< 1 . 


(BIO) 


equivalent to F(Ef,p) > 0, where we have defined 
F(E f ,p) = 4p 2 (E 2 - m 2 ) 

- [El - mix -p 2 + m 2 - E 2 f ] 2 . 


(Bll) 


For a given value of the absolute momentum the roots of 
F(- ,p) read for the plus-sign case 0 = Ef + E b + po , and 
therefore El = ( Ef +p 0 ) 2 : 


£4 ip,v n ) = g 


irqq 


[ d \ 

1 

E - 1 [ 

-^maXjinin — 2?77^ 

1 (2tt)3 

[2 E b 



E\ 


Z 2 


(■ E 2 -E 2 + u 2 n )Z 3 


Zl - + - 

vl E 2 _+vl) 2 w n {E\ + vl)(E 2 _ + i/2) 

(B 6 ) 

The non-vanishing imaginary parts of £j are induced by 
their pole structure: 


( to m — 2 m 2 ) \Jm 2 + p 2 


±P m M \Jm-M - 4 m 2 


(B12) 


The range of integration, Ef £ [E m i n , E max ], depends 
therefore linearly on the external quark momentum: 


■Ehiax - -Emin = P ^ \/m^ - 4m 2 . 


lim Im 

£—>•0 




(B13) 


= ZttS(x 2 — Pq) 

v n \-+—ipo+e 

IT 7 

= if— ( s ( x - Po) + S(x + Po)) ■ 

2po 


(B7) 


In the limit of a vanishing external quark momentum the 
range of integration collapses to one single point: 


E, 


m 


This means for the Z\ term: Ef + E b ± po = 0, where 
only the minus sign can be realized. For the Z 2 term, 
Ef — E b ± po = 0, both signs can be realized for the 
time being. We will see that only the plus-sign case 
contributes to the (on-shell) imaginary parts, so there 
is just one contribution from Z 2 . Later, the Z 3 term is 
considered separately. We start with the first two terms 
Z £ {Zi, Z 2 }. Using the identify (B71 we find the follow¬ 
ing structure when evaluating £ 0i 3 (p, v n ) from Eq. (B2) 
after analytical continuation has been carried out: 


max. min 


p—0 


_M 

2 m 


— m > m . 


(B14) 


We emphasize that the whole discussion is only valid for 
temperatures above the Mott temperatures Tm, where 
the pion mass is at least twice the constituent-quark 
mass. This constraint can be seen explicitly from the 
integral boundaries Eq. (B12 1 . We have already intro- 


aq nZ 1 

(2tt) 3 2po 2 E b E f 
f d 3 g nZ 


5(E b - (*)) 


(*,W <E *- ( * )2) 

1 /*oo j_2 


= 2tt 


= 2tt 


I dC / 

-1 Jo 

r dE f 


dg q 2 nZ 
(2n ) 3 2 p 0 E f 
irZ 


(B 8 ) 


«(0-(*r) 


(27r) 3 4p 0 p 


©a-n, 


where £ = cosfh In order to carry out the integral over 
the delta function we have used 


E b ~ m M + 
dE'l 


dt 


(p - q ) 2 = Tom +p 2 +q 2 - 2 pq£ , 
= 2 pq , 


duced the Mott temperature when discussing thermal 
quark and meson masses, where Tm ~ 212 MeV have 
been found in the case of vanishing quark chemical po¬ 
tential. Note that this discussion remains valid also in 
the chiral limit, where the current-quark mass is set to 
zero, mo = 0. In this case, the pion mass vanishes in 
the Nambu-Goldstone phase at low temperatures but it 
is finite when chiral symmetry is restored for large tem¬ 
peratures. 

We conclude that under the condition mM > 2m, i.e. 
for T > Tm, the phase space is always non-empty and 
compact: 0 ^ [E m - ln , E max ] C [m, 00 ). This fact implies 
that the shear viscosity 77 will evaluate to some finite 
result in this temperature region. However, we have also 
derived the following substitution rule 


J (-)e(i-e 2 ) =/_“ ffi(-), (B15) 


(27r) 


(B9) 


which leads finally to a well-conditioned one-dimensional 
numerical integral. 
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For the sake of completeness, we also mention the 
minus-sign case, i.e. 0 = E b + Ef — pp . If we plug in 
E 2 = (Ef — p 0 ) 2 into the condition (BIO I then the phase 
space simply vanishes for any incoming quark momen¬ 
tum, since the range of integration would be restricted 
to negative energies in the fermion loop: 


T^min — i -^max — ^min 


(B16) 


We can therefore conclude that only the plus-sign case, 
E b = Ej +p 0 , allows for an on-shell condition for the 
mesonic fluctuation. Knowing this we can now continue 
with the third term, Z 3 , in Eq. (B2): 


i v n Z 3 


\v n >-4po+i£ 


Im [(Ef + E b y + ul\ ((Ef ~E b y + ul\ 

= Po nZ 3 5 (((Ef + E b ) 2 - pl\((Ef - E b ) 2 - - 2 
P 0 ttZ 3 


2 po 

nZ 3 

AE f 


6 (((Ef + E b ) 2 -p 2 0 )(E b -Ef-p 0 }) 

"-V-' 

=4 Ef Eb 

6 (E 2 -(E f +p 0 ) 2 ). 


(B17) 

Note that due to the \v n factor in the first line, the po 
terms cancel in the final result. As done in the calculation 


(B 81 the momentum integral can be performed: 


aq nZ 3 
(2tt ) 3 4Ef 


5(E 2 - (Ef + p 0 ) 2 ) 


= 2 7T 


d Ef ttZs 


(B18) 


e(i-a. 


(27t) 3 8 p 

Combining all contributions, we find with J -4 = — - 


Im E 0 i 3 , 4 (p,- ipo) 

q2 /* -E'max 

= -ff L / dE f Eo, 3 A [n B (E b )+n^(Ef)\, 

E J i? m in 

(B19) 


as it has been stated in Eq. (391. 


Appendix C: Details of the off-shell calculation 


Due to the symmetry properties of the quark spectral 
function, p(—po,p,T,p) = —p(po,p,T,—p), we can sim¬ 
plify our discussion and restrict the off-shell energy to 
non-negative values po > 0 . 


We return to H 0 ,3 given in Eq. (B2) and decompose 
into partial fractions for convenience: 


S 0,3 (p, -ipo) = flMqq J 


d 3( Z ^o,3 


(27r) 3 AE b E f 

1 - rip(Ef) + n B (E b ) n B (E b ) + rip(Ef) 
Ef + E b + po Ef — E b + po + is 

n B (E b ) + rip (Ef) 1 + n B (E b ) - n^(E f ) 
Ef - E b - po - ie Ef + E b - p 0 - ie 


(Cl) 


As in the on-shell discussion, taking its imaginary part 
probes the pole position of the partial fractions intro¬ 
ducing four cases ±E b = Ef ±p 0 . The fraction in the 
first line introduces E b = —Ef — po < 0 which can be ex¬ 
cluded immediately. The remaining three cases are de¬ 
noted as follows: 


Case I: E b = Ef+po, 

Case II: E b = Ef — p 0 , (C2) 

Case III: E b = p 0 — Ef . 


As mentioned we restrict the discussion to po > 0. Car¬ 
rying out the d 3 g integral introduces again the restriction 
|£| < 1 with £ = cos 8 denoting the polar angle: 


_ 1 < E b (p,p 0 ) 2 — TO^ — p 2 — q 2 < 
2 pq 


(C3) 


equivalent to F(Ef,p,po) > 0, where we have defined 
E(Ef,p) = Ap 2 (E 2 — m 2 ) 


- [E b (p,p 0 ) 2 - Wm - P 2 + m 2 - 


(C4) 


In the off-shell case p > 0 and po > to are independent 
of each other. In the following we evaluate the three- 
dimensional integral (Cl) ensuring |£| < 1 by applying 


the three cases for the relation between quark and meson 
energy. 

Case I. This is the only case that can be realized on- 
shell: E b = Ef+po- The following two conditions have to 
be fulfilled: (i) Ef > m and (ii) Ef > tom~Po, which can 
be summarized in Ef > max(m,mM — Po)- Evaluating 
the condition £ < 1 we find 


E(E f ,p,po) > 0 <:=> ~4s(E f -E-)(E f - E + ) > 0 , 

(C5) 

where we have introduced s = p 3 — p 2 and 

£ Po , ( m M ~ m2 )Po 

2 2s (C6) 

± — V[s - (m + to m ) 2 ] [s - (to - to m ) 2 ] ■ 


These roots of F( Ef, •) are generalizations of E, 
introduced in Eq. (|B12 1 . One finds indeed 




= E, 


s=Pq— p 2 —m 2 >0 


max,min 


max,mm 


(C7) 


Note that in contrast to m < E m - ln < E m ax , the off-shell 
roots are not ordered that simply. Dependent on s > 0 
or s < 0 one has E_ < E + or E + < E_, respectively. 
In addition, it might happen that one or even both roots 
are negative as we will see. 

First, we consider the case s < 0 which leads to a 
convex-up parabola F(Ef ,-,-) with possible integration 
ranges Ef < E + and Ef > E_. In general one has 
to distinguish additionally the two cases m < tom and 
to > tom, but right now we find for both cases 


E + < —Tom — Po < 0 , E_ > max(m, tom- P o ) • (C 8 ) 
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(a) s < 0 



(b) s > 0 with m < my[ and po < yj(m — m-M) 2 + P 2 

Figure 9. Summary of integration ranges (gray boxes) for the 
off-shell imaginary parts of £o, 3 , 4 (p, —ipo) for Case I. On-shell 
only the case (b) can be realized. 


For Case I with s < 0 we have the range of integration 
Ef > E_ as sketched in Fig. [ 9 J) a). 

Now consider the case s > 0 with a concave-down 
parabola F(Ef,-,-). The possible integration range is 
E_ < Ef < E + . This time, the roots are not automati¬ 
cally real numbers, but for (m — tom) 2 < s < (to + tom) 2 
they become purely imaginary and have to be excluded. 
The first option p 0 < yj(m — Tom) 2 + p 2 leads to 

to < tom : E + > E- > ma x(to, tom — Po) , 

~ ~ (C9) 

to. > tom : E- < E + < —to — po < 0 . 

Therefore, the case to > tom cannot be realized and only 
for to < tom the full range of integration is accessible. 
We summarize this case in Fi g.^b). Having s > 0 there 
is the second option po > \J (m + tom) 2 + V 2 for which 
one has for both cases to < tom and to > Tom: 

E + < —m , E_ > m m — Po ■ (CIO) 

We conclude Tom — Po < E_ < E + < — m < 0, hence 
this option is excluded and the discussion of Case I is 
completed. 

Case II. Evaluating the condition |£| < 1 using Et = 
Ef — po leads to 

F{E f ,p,po) > 0 —4 s(E f +E-)(E f +E+) > 0 , (Cll) 



(a) s < 0 



(b) s > 0 with m > mM and po < \/(m — mm) 2 + P 2 

Figure 10. Summary of integration ranges (gray boxes) for 
the off-shell imaginary parts of Eo, 3 , 4 (p, —ipo) for Case II. 


hence —E± are the roots of F(Ef ,-,•)• We follow the 
same systematic path as before: 

Consider first the case s < 0 implying again a convex- 
up parabola with E_ > E + . This means —< — E + , 
providing two possible integration ranges Ef < — E_ and 
Ef > —E + . One finds: 

to < to m : — E + > to m + Po > 

- E- < -to m + Po < rn M + Po , 

to > tom : — E + > tom + Po , 

— < —m + po < -tom + Po < % + Po ■ 

(C12) 

In conclusion we find the range of integration as shown 
in Fig. [l0];a), again without any restriction on the quark 
and meson masses. 

Now consider s > 0. This time the possible range of 
integration is — E+ < Ef < — E-. For the option po < 
\J (to — tom) 2 + p 2 we find 

to < tom : - E + < —E_ < —tom + Po < % + Po > 

to > Tom : — E_ > —E + > tom + Po • 

(C13) 

Using the constraint Ef > m m + Po the case m < tom 
is excluded and only m > tom is possible. The option 
Po > \J (to T Tom) 2 + P 2 , for both cases to < Tom and 
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(a) s > 0 with p o > y/{m + m M ) 2 + p 2 


with 


J 1 


0(p-Po) 


H - po + T In 


M^~) 

«b(-E- +p 0 ) 


+ #(Po - p)0(m M - m)6{yj(m - m M ) 2 + P 2 - Po) 

x T in n F^-^ nB ^+ + p °) 
np(-E + )n B (-E- + Po) 

(C17) 


Figure 11. Summary of integration ranges (gray boxes) for 
the off-shell imaginary parts of Eo, 3 , 4 (p, —ipo) for Case III. 


to > Tom, leads to: 


- E + < -E_ < -to m + Po < m M + Po , (C14) 

which excludes this case because Ef > max(m, Tom +Po) 
must be ensured. This case is illustrated in Fig. Un¬ 
case III. The final case, E b = po — Ef, leads to the 
two conditions (i) Ef > to and (ii) Ef < p 0 — Tom- From 
this we get po > to + tom- It is important to realize that 
Eb in this case is just the negative of the condition used 
in Case I. Therefore, all contributions present for Case I 
cannot be realized for Case III. It remains to check the 
case s > 0 in combination with po > yj (to + tom) 2 +P 2 - 
We have — E + < Ef < —E_ as possible integration range 
and find (cf. the related discussion for Case I): 

- E + > to , -E_ < po - to m , (C15) 


J =e(p-p 0 ) 


H - po + T In 


n B (-E+ - po) 

n F (~E+) 


+ 6 {Po ~ p)0(m - to m )6'(v / (to - to m ) 2 +p 2 - Po) 
x T In n F(~E-) u b(-E+ ~ Po) 


~{-E + )n B (-E _ -p 0 ) 


(C18) 


J m = 9(p 0 - yj(vn + m M ) 2 + p 2 ) 

x T In n F(E+)‘>ME- +Po) 
n.p{E_) n B (E + +p 0 ) 


(C19) 


Note that for J lu the condition 9(p 0 —p) just follows from 
9(po — yj (to + tom) 2 + p 2 ), therefore this 0-function can 
be omitted. 


which is valid for both to < Tom and to > tom- In con¬ 
clusion, there is only one contribution to the imaginary 
part for Case III as shown in Fig. m 

Combining now all three cases, the off-shell imaginary 
part of So(p, —ipo) can be calculated immediately. The 
rather lengthly result read^] 


Im E ° ff = fS {l dEf Hb(£/+po) - M^/)] 

+ / dE f [n B (Ef -Po) + r4(Ef)] 

J ii 

+ d E f [l + n B (p 0 - E f ) - n+ (E f )] | 

( J 1 + J 11 + J HI ), 


0Mqq ( T \ j 7 H , 7 IIINj 

167Tp 


(C16) 


Next we present the off-shell result for the imaginary 
part of S 4 (p, —ipo), performing again a partial-fraction 
decomposition of Eq. (B 61 : 


S 4 (p, ipo) <?Mqq J ( 27r)3 ApoEbEf 

x (1 ^ rip(Ef) +n B (E b ))(E b +p 0 ) 

L E f + E b + po 

{n B (E b ) + np (Ef))(p 0 - E b ) 

Ef — E b + po + ie 
(n B (Eb) +n£(Ef))(p 0 - E b ) 

Ef - E b - po - ie 
| (1 + n B (E b ) - n^{Ef))(p 0 - E b ) 

Ef + E b - po - ie 


(C20) 


‘‘ The minus signs for Case I is due to the pole description+ie 
instead of —ie for Case II and Case III. For E b we have always 
inserted the corresponding relations to Ef and po as defined in 
Eq. JC2l. 


In comparison to Eq. (B21 there is the factor po in the 


denominator and also combinations of E b and p 0 in the 
numerators, but X 4 features the very same pole structure 
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as discussed before. Therefore we find immediately: 


Im E f = 


9m 


qq 


167rppo 


X {Ii dEf ^ Ef ^~ nB ^ Ef + p °) “ n F( E f)] 

+ f d E f E f [n B (E f -p 0 ) + n£(E f j] 

J ii 

+ dE f E f [l + n B {po - E f ) - (E f )] j 


9m 


qq 


(. K l + K n + K m ) . 


16npp 0 

Introducing the two auxiliary functions 


(C21) 


g ± (E) = T 


2 


we find: 


y + Li 2 (l - + U 2 

(C22) 


K U = 0(p~Po) |^(M 2 ~Po) +E+T In n%(-E+) 

+p 0 T In n B (-E + - p 0 ) + CT (-£?+)} 

+ 0{po - p)9(m - m M ) 8 (\J{m - m M ) 2 + P 2 - Po) 

x L 0 T In M-E+-PO) + fg- { _E + ) _ g-{-E_)\ 

{ n B (-E^-p 0 ) v ' 

+E+T In n+(-E+) - E_T lnn+(-£_)} , 

(C24) 

K m = 0(p o - yj (m + to m ) 2 +P 2 ) 

x {- p 0 T In n B(^++Po) _ (g+(E_) - g+(E + )) 

[ n B (E_+p 0 ) v > 

+E_T \nrip(E_) — E + T lnnp (£’+)| . 

(C25) 

Note that the expression for I \ 111 is just the negative of 
the second contribution to K 1 . 

Having derived results for Im Eg ff and Im E^®, the re¬ 
maining integral for E 3 can be performed easily, since 
all building blocks have been prepared. The main obser¬ 
vation is that E 3 splits into two parts: the first one is 
independent of Ef, the second one introduces the same 
Ef dependence present in the calculation for E 4 : 

2p 2 E 3 = m 2 M +p 2 +q 2 -E 2 b 

{ m M + P 2 ~ Po ~ 77,2 ~ 2Efp 0 for Case I , 
m M +P 2 ~ Pq — m 2 + 2Efp 0 for Case II, III . 

(C26) 

We find therefore 


K 1 = 0{p - Po) | ^{p 2 3 - Po) - E-T lnn F (£_) 

-PoT In n B {E- +p 0 ) + f? + (-E_) j 

+ 9{p 0 - p)0(rriM - m)0{yj (m - m M ) 2 + P 2 - Po) 

x Lt hr ME++P0) + (g+{E_) - g+(E + )) 

[ n B (E_+p 0 ) v 1 

+E+T lnn F {E + ) — EJT lnn F (£_)j , 

(C23) 


Im Ef = ^ [-"b (Ef +Po) - np(Ef)] 

+ I d Ef El 1 ' 111 [n B {E f - po) + n+(E f )\ 

J 11 

+ d E f El 1 ’ 111 [1 + n B {po - E f ) - n+(E f )] J . 


(C27) 

Inspecting the definitions for J x in Eq. (C16l and A' x in 
Eq. (C21), we find the relation: 


Im Ef = im S f + P 1 Im Ef 

2 p z p- 


(C28) 
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